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FIGURES 


Frequency  spectra  of  spread-spectrun  signals  . 

Effects  of  Pseudonoise  Despreading  on  noise  and  jamming 


1 .  INTRODUCTION 

Success  in  modern  military  warfare  depends  heavily  on  maintaining 
adequate  and  secure  lines  of  communication.  Thus,  the  attempted 
interception  and  jamming  of'  communications  by  the  enemy  seems 
inevitable,  and  measures  must  be  taken  by  the  communicators  to  ensure 
that  their  lines  of  communication  remain  open  and  secure.  One  method  of 
decreasing  the  vulnerability  of  communication  systems  to  interception 
and  jamming  is  to  employ  spread-spectrum  techniques.  These  techniques 
provide  resistance  to  jamming  and  interception  by  distributing  the 
transmitted  energy  over  large  bandwidths.  The  energy  density  of  the 
transmitted  waveform  is  therefore  decreased,  yielding  a  low  probability 
of  intercept  (LPI)  system.  In  addition,  if  unsophisticated  jamming  of 
at  least  a  significant  portion  of  the  signal  bandwidth  is  desired,  then 
the  bandwidth  over  which  a  jammer  must  operate  is  increased,  forcing  the 
jammer  to  dilute  its  available  power.  This  dilution  results  in 
decreased  jamming  energy  density  and,  consequently,  decreased  jamming 
effectiveness. 

The  quantification  of  the  decrease  in  jamming  effectiveness  and  the 
corresponding  increase  in  the  intelligibility  of  the  communication 
signal  is  not  a  simple  procedure,  however.  In  fact,  signal 
intelligibility  is  a  highly  complex  function  of  the  signal  and  jamming 
parameters,  atmospheric  and  terrain  conditions,  and  a  number  of  other 
factors.  Voice  communication  intelligibility  is  particularly  difficult 
to  theoretically  quantify  because  it  is  heavily  dependent  on  the 
communicators'  speaking  and  listening  abilities  as  well  as  on  message 
content.  Digital  data  intelligibility,  on  the  other  hand,  is  somewhat 
simpler  to  quantify  in  terms  of  bit  error  rates,  but  the  use  of  bit¬ 
interleaving  techniques  or  error-correcting  codes  in  the  communication 
system  can  complicate  the  determination  of  word  error  rates.  Spread- 
spectrum  techniques,  because  they  tend  to  decrease  transmission  errors 
in  the  presence  of  jamming  or  other  interference,  further  complicate  the 
determination  of  digital  message  intelligibility. 

This  paper  describes  a  computer  program  which  analyzes  the 
performance  of  digital  communication  systems  that  employ  certain  types 
of  spread-spectrun  techniques  in  the  presence  of  interfering  signals. 
The  program  computes  the  bit  or  word  error  rate  as  a  function  of  the 
communication  signal  parameters  ( spectrvm  spreading  technique,  error- 
correcting  codes,  data  modulation,  bandwidth,  etc.)  and  the  interference 
parameters  (thermal  noise,  jamming  power,  jamming  modulation,  jamming 
bandwidth,  etc.).  The  program  output  consists  mainly  of  plots  of  bit  or 
word  error  rates  as  functions  of  one  or  more  of  the  signal  or 
interference  parameters,  as  determined  by  the  user.  The  plots  are 
produced  interactively  on  a  Tektronix-type  storage-tube  graphics 
terminal  connected  to  an  IBM  370/168  computer  which  executes  the 
program.  The  use  of  interactive  computer  graphics  provides  rapidly 
produced  results  and  high  user-computer  interaction. 
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The  computer  model,  as  well  as  the  theoretical  treatment1  on  which 
the  model  is  based,  assumes  that  the  spread-spectrum  receiver  has 
acquired  the  intended  signal  and  that  the  transmitter  and  receiver  are 
in  perfect  synchronization.  Thus,  the  effects  of  jamming  on  signal 
acquisition  and  synchronization  are  not  modeled  in  this  program. 

Since  the  computer  model  described  in  this  report  is  based  on  the 
theoretical  treatments  of  Torrieri1  '2  and  since  this  report  frequently 
references  those  treatments,  users  of  this  report  should  have  CH/CCM-78- 
2  and  CM/CCM-79-9  readily  available. 


2.  PROGRAM  OVERVIEW 

The  computer  program  described  in  this  paper  was  designed  to  analyze 
and  compare  digital  communication  systems  employing  two  types  of  spread- 
spectrum  techniques:  frequency  hopping  (FH)  and  pseudonoise  (PN,  also 
known  as  direct  sequence  or  DS),  in  environments  where  considerable 

interference  is  present. 

Frequency  hopping  is  defined  as  the  periodic  changing  of  the 

communication  signal-carrier  frequency  throughout  a  wide  band  of 
frequency  choices.  A  system  that  changes  frequency  at  a  rate  higher 
them  the  data  rate  of  the  digital  information  to  be  transmitted  is 
called  a  fast-frequency  hopping  system.  A  system  that  hops  slower  than 
the  information  rate  is  called  a  slow- frequency  hopping  system. 

Although  the  bandwidth  of  a  frequency  hopping  signal  at  any  given 
instant  during  transmission  is  small  ( aproximately  equal  to  the  data 
rate  in  slow-hopping  systems  or  the  hopping  rate  in  fast-hopping 

systems) ,  the  total  bandwidth  occupied  by  the  signal  over  a  long  period 
of  time  is  much  greater  than  this  "instantaneous"  bandwidth  because  of 
the  hopping  carrier  frequency  (see  fig.  1).  The  resistance  of  frequency 
hopping  signals  to  intercept  and  jamming  derives  from  the  large  total 
bandwidth  generated  by  the  hopping  systems.  In  order  to  be  effective 
against  most  practical  hopping  systems,  a  potential  interceptor  or 
jammer  must  either  follow  the  hopping  signal  in  frequency  or  increase 
its  bandwidth  to  encompass  at  least  a  significant  portion  of  the  hopping 
bandwidth.  The  former  tactic  is  a  difficult  task,  especially  for 
systems  with  high  hopping  rates;  the  latter  tactic  results  in  a 
reduction  in  the  interceptor's  signal-to-noise  ratio  (making 
interception  more  difficult)  and  also  in  the  amount  of  jamming  power 
present  within  the  relatively  narrow  "instantaneous"  bandwidth  of  the 
hopping  signal . 

1 D.  J.  Torrieri,  Frequency  Hopping  in  a  Jamming  Environment,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
( December  1978 ) . 

2D.  J.  Torrieri ,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  DARCOM  Report 
CM/CCM-79-9  (December  1979). 
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Figure  1.  Frequency  spectra  of  spread  spectrum  signals. 


Pseudonoise  modulation  is  defined  in  this  report  as  the  periodic 
alteration  of  the  phase  of  the  transmitted  waveform  achieved  by  phase- 
shift-keying  the  carrier  after  data  modulation  has  been  completed.*  If 
the  frequency  at  which  phase  alterations  occur  is  considerably  higher 
than  the  transmitted  data  rate,  then  the  null-to-null  bandwidth  of  the 
transmitted  signal  is  increased  to  approximately  twice  the  frequency  at 
which  phase  alterations  occur  (see  fig.  1).  The  resistance  of 
pseudonoise  signals  to  intercept  and  jamming  derives  from  the  large 
bandwidth  generated  by  the  pseudonoise  system.  Since  the  transmitted 
energy  is  spread  over  a  large  bandwidth,  the  spectral  power  density  of 
the  signal  is  reduced,  and  in  some  systems  is  below  the  noise  level. 
The  pseudonoise  signal  is  noise-like,  so  a  potential  interceptor  would 
have  great  difficulty  even  detecting  the  presence  of  the  signal.  For  a 
given  amount  of  total  jamming  power  in  the  pseudonoise  bandwidth,  the 
resistance  of  the  pseudonoise  signal  to  jamming  does  not  depend  greatly 
on  the  bandwidth  of  the  jamming  because  of  the  despreading  function  in 

*  Although  other  forms  of  pseudonoise  modulation  exist,  phase-shift 
keyed  pseudonoise  modulation  is  the  only  type  cons idered  in  this  report 
because  it  is  the  most  commonly  used  type  of  pseudonoise  modulation . 
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the  receiver.  This  despreading  function  merely  demodulates  the  phase- 
shift-keyed  pseudonoise  modulation  imposed  on  the  signal  by  the 
transmitter ,  thereby  restoring  the  signal  to  its  original  form. 
However,  the  despreading  function  is  actually  the  inverse  of  the 
spreading  function,  and  the  application  of  the  despreading  function  to 
any  signals  other  than  the  intended  pseudonoise  signal  will  result  in 
the  spreading  of  those  signals.  Thus,  narrowband  jamming  is  effectively 
spread  by  the  despreading  function,  resulting  in  the  wideband  noise- like 
spectrun  in  figure  2.  wideband  noise  jamming  is  not  affected  by 
despreading  because  of  the  random  nature  of  this  type  of  jamming  (that 
is,  further  randomization  of  an  already  random  process  has  no  effect) . 
Thus,  wideband  jamming  retains  its  reduced  effectiveness  (caused  by  its 
wide  bandwidth  and  reduced  power  density) ,  while  narrowband  jamming  is 
altered  by  the  despreading  function  to  look  like  wideband  jamming. 
Additional  information  about  frequency  hopping  and  pseudonoise 
modulation  is  available .  ^ 

In  both  frequency  hopping  and  pseudonoise  communication  systems, 
there  must  be  some  method  which  tells  the  transmitter  what  frequency  to 
change  to  or  when  to  change  phase.  This  method  must  also  be  known  to 
the  receiver  so  that  it  can  follow  the  transmitted  signal  in  frequency 
or  phase  in  order  to  derive  the  transmitted  data  as  they  were  before 
spectrum  spreading  was  imposed.  The  method  used  in  most  frequency 
hopping  and  pseudonoise  systems  uses  a  sequence  of  binary  digits  called 
a  spreading  code.  The  spreading  code  is  normally  (although  not 

necessarily)  a  pseudo-random  generated  code,  so  that  interception  of  a 
portion  of  the  communication  signal  will  not  allow  the  determination  of 
the  code-generation  scheme.  In  typical  frequency  hopping  systems,  the 
spreading  code  bit  generated  for  each  hopping  period  is  shifted  into  a 
shift  register  containing  a  number  of  previously  generated  spreading 
code  bits,  and  the  contents  of  the  shift  register  are  used  to  determine 
the  next  frequency  in  the  hopping  sequence.  In  typical  binary 
pseudonoise  systems,  the  new  spreading  code  bit  is  compared  to  the 
previous  bit,  with  a  change  in  logical  state  between  these  two  adjacent 
bits  causing  phase  alteration  of  the  carrier  signal.  The  rate  at  which 


*£>.  J.  Torrieri ,  Frequency  Hopping  in  a  Jamming  Environment ,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
(December  1978). 

2D.  J.  Torrieri ,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  Report  CM/CCM-79-9 
(December  1979). 

3 Spread  Spectrum  Operator's  Handbook,  National  Security  Agency,  NSA 
Report  W32-228-78  (30  May  1978). 

4R.  C.  Dixon,  Spread  Spectrum  Systems,  John  Wiley  and  Sons,  Inc.,  New 
York  (1976). 
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the  spreading  code  is  generated  determines  the  hopping  rate  in  a 
frequency  hopping  system,  but  in  general  does  not  determine  the  total 
bandwidth  of  the  hopping  signal.  (Although  the  hopping  rate  does  have 
an  effect  on  the  "instantaneous"  signal  bandwidth  and  thus  may  alter  the 
total  bandwidth,  this  alteration  is  not  significant,  except  in  systems 
with  very  high  hopping  rates  and  relatively  narrow  total  bandwidths.) 
The  total  bandwidth  of  a  frequency  hopping  system  is  determined 
primarily  by  the  range  of  possible  frequency  choices  (or  "channels" } 
allowed  by  the  system.  As  previously  mentioned,  the  null-to-null 
bandwidth  of  a  pseudonoise  system  depends  on  the  spreading  code  and, 
specifically,  the  bandwidth  is  approximately  twice  the  spreading-code 
generation  rate. 


(A)  BEFORE  DESTREADMC 


(B)  AFTER  DE8PREADM6 


Figure  2.  Effects  of  pseudonoise  despreading  on  noise  and  jamming. 
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It  is  intuitively  obvious  that  the  use  of  frequency  hopping  or 
pseudonoise  modulation  on  a  signal  decreases  the  vulnerability  of  the 
signal  to  narrowband  interference  by  reducing  the  signal's  power  density 
and  increasing  its  bandwidth.  However,  it  is  difficult  for  on  bo 
intuitively  determine  exactly  how  much  improvement  either  of  these 
techniques  might  offer  in  any  particular  communication  application. 
Theoretical  treatments  of  the  effects  of  frequency  hopping  and 
pseudono.se  modulation  have  been  published.1'2  The  treatment  uses 
strictly  analytical  models  of  these  spread-spectrum  techniques  to  derive 
expressions  for  bit  and  word  error  rates  in  terms  of  the  characteristics 
of  the  spectrum  spreading  techniques  being  used,  the  information  being 
transmitted,  and  the  interference  that  is  present. 

The  program  described  in  this  paper  is  a  computer  adaptation  of 
portions  of  the  theory  presented  by  Torrieri.1,2  The  program  models  the 
following  communication  and  interference  characteristics. 

( a)  Spread- spec trun  techniques: 

1.  Fast-frequency  hopping — models  equation  (14)  of  CM/CCM- 

78-2. 1 

2.  Slow-frequency  hopping — models  equations  (20)  and  (59) 

of  CM/CCM-78-2. 

3.  Pseudonoise  (direct  sequence) — models  equations  (54)  and 

(55)  of  CM/CCM-79-9.2 


( b)  Data  modulation  techniques  (frequency  hopping  only): 

1.  Binary  frequency  shift  keyed  (binary  FSK) — This  type  of 
modulation  uses  two  different  frequency  channels  for  each  transmitted 
information  bit;  transmission  at  one  frequency  denotes  a  logical  "1," 
while  transmission  at  the  other  frequency  denotes  a  "0."  Demodulation 
of  the  received  signal  is  done  through  a  comparison  of  the  energy  in  the 
two  channels;  the  higher  energy  is  selected.  Thus,  the  receiver  does 
not  require  any  threshold  circuitry.  The  advantages  of  this  type  of 
modulation  are  that  it  requires  relatively  simple  and  inexpensive 
hardware,  and  it  does  not  require  phase  coherence  of  the  carrier  from 
bit  to  bit.  Its  main  disadvantage  is  that  it  requires  two  channels  for 
each  transmitted  bit  and  therefore  uses  twice  the  bandwidth  of  other 
modulation  techniques. 


1D.  J.  Torrieri ,  Frequency  Hopping  in  a  Jamming  Environment ,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
(December  1978). 

2D.  J.  Torrieri,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  Report  CM/CCM-79-9 
( December  1979 ) . 


2.  Coherent  phase  shift  keyed  (CPSK).  This  type  of  modulation 

encodes  the  information  bit  to  be  transmitted  in  the  phase  of  the 

carrier.  Hie  receiver  demodulates  this  signal  by  measuring  the  received 
phase  and  comparing  it  to  a  reference;  an  in-phase  or  out-of-phase 
condition  determines  the  value  of  the  received  bit.  Hie  advantages  of 

this  type  of  modulation  are  its  reduced  bandwidth  requirements  (half 

that  of  binary  FSK)  and  its  higher  resistance  (with  respect  to  FSK)  to 
bit  errors  due  to  interference.  Its  main  disadvantages  are  that  it 
requires  complex  circuitry  in  the  receiver  for  thresholding  and  phase 
measurement  and  that  carrier  phase  coherence  must  be  maintained  by  the 
transmitter  when  generating  the  information  bits . 


(c)  Error-correcting  codes:  The  program  allows  the  use  of 

repetition  coding,  where  each  information  bit,  is  transmitted  more  than 
once  and  at  more  them  one  hopping  frequency.  The  receiver  uses  majority 
rule  to  determine  the  correct  value  of  the  information  bit.  Hie  program 
also  allows  Hamming-type  error-correcting  codes,  such  as  (7,4)  coding. 
Hiis  type  of  coding  converts  a  4-bit  information  word  to  a  7-bit 

transmitted  bit  stream.  The  receiver  is  able  to  detect  up  to  two  bit 
errors  in  the  7-bit  stream,  and  can  correct  one  bit  error  in  the  stream, 
thereby  correctly  restoring  the  4-bit  information  word  in  spite  of  the 
single  bit  error.  Hie  program  accepts  any  type  of  (m,n)  coding;  the 
user  need  supply  only  m  (the  total  number  of  bits  after  coding),  n  (the 
number  of  information  bits  per  word  before  coding) ,  and  the  maximum 

number  of  bit  errors  that  the  code  can  correct  in  each  m-bit  block. 

(For  [7,4]  coding,  these  numbers  are  7,4,  and  1,  respectively.) 

( d)  Jamming  techniques; 

1.  Wideband  Gaussian  noise  jamming  with  variable  bandwidth. 

2.  Narrowband  (continuous-tone)  jamming. 

3.  Narrowband  repeater  jamming  (sometimes  called  follower 

jamming) . 

4.  White  noise  repeater  jamming. 

The  user  of  the  program  also  specifies  the  following  parameters. 


(e)  Information  bandwidth:  The  bandwidth  of  the  information  signal 
before  spread-spectrun  techniques  are  applied.  This  value  is  also  used 
by  the  program  as  the  frequency  hopping  channel  bandwidth. 

(f)  Total  bandwidth:  The  bandwidth  of  the  transmitted  signal  after 
spread-spectrum  techniques  are  applied. 

( g)  Number  of  hopping  channels:  The  total  number  of  frequency 

choices  (or  pairs  of  frequency  choices  for  binary  FSK  modulation) 
available  for  frequency  hopping.  The  number  may  either  be  specified  by 
the  user  or  derived  in  the  program  by  dividing  the  total  bandwidth  (f. 
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above)  by  the  information  bandwidth  (e,  above) .  The  latter  derivation 
results  in  the  use  of  the  maxim  van  number  of  hopping  channels  available 
in  the  total  bandwidth. 


(h)  Signal -to-noise  ratio  (SNR):  The  ratio  in  decibels  of  the 

total  signal  power  to  the  thermal  noise  power. 

(i)  Jamming-to-signal  ratio  (JSR):  The  ratio  in  decibels  of  the 

total  jamming  power  present  over  the  total  spread-spectrun  bandwidth  to 
the  total  signal  power .  In  other  words,  JSR  is  the  jamming-to-signal 
ratio  that  would  exist  if  no  spec trim  spreading  techniques  were  used  and 
if  the  total  jamming  power  were  concentrated,  within  the  information 
bandwidth . 

The  program  does  not  consider  signal  acquisition  and  synchronization 
in  its  analyses  of  spread-spectrum  communication  systems. 
Synchronization  is  the  procedure  that  aligns  the  pseudo- random  codes 
used  in  generating  the  transmitted  spread- spectrum  signal  with  the  codes 
used  in  the  despreading  portion  of  the  receiver  to  allow  proper 
despreading  to  occur .  The  problem  of  synchronization  in  spread-spectrum 
systems  is  similar  to  the  problem  of  synchronization  in  enciphered 
systems  because  enciphered  systems  also  employ  pseudo- random  codes  that 
must  be  properly  aligned  in  the  receiver  before  deciphering  can  be 
accomplished.  Since  transmitter-receiver  synchronization  is  required 
before  any  communication  may  take  place,  this  is  often  considered  one  of 
the  more  vulnerable  links  in  a  system.  In  particular,  hostile  jamming 
may  be  targeted  on  a  communication  system  specifically  to  deny  or 
disrupt  the  synchronization  process.  Although  this  process  is  an 
important  part  of  any  spread-spectrum  system,  it  is  difficult  and 
complicated  to  analyze  or  model  because  of  the  many  variables 
involved.  Therefore,  both  the  program  described  in  this  paper  and  the 
theoretical  treatments  of  Ttorrieri^'^  on  which  this  program  is  based 
assume  that  transmitter-receiver  synchronization  has  been  accomplished 
and  is  maintained  throughout  the  duration  of  the  communication  in 
question. 


3.  PROGRAM  DESCRIPTION 

The  computer  program  described  in  this  report  employs  interactive 
computer  graphics  to  provide  improved  user-computer  interaction  and  a 
simple,  easily  understood  graphical  display  of  the  results  of  an 
analysis.  Dialogue  between  the  user  and  the  program  is  generally  in  a 
question-answer  format.  The  user  is  first  supplied  with  a  short 

1D.  J.  Torrieri,  Frequency  Hopping  in  a  Jamming  Environment,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
(December  1978). 

2-D.  J.  Torrieri ,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  Report  CM/CCM-79-9 
(December  1979) 
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overview  of  the  program.  The  user  is  then  asked  a  series  of  questions 
concerning  the  characteristics  of  the  ccomunication  system  and  the 
interference  environment  to  be  analyzed.  The  program  displays  each 
question,  and  then  waits  after  each  question  until  the  user  has  provided 
an  answer.  Some  questions  require  only  a  yes  or  no  answer;  others 
require  the  user  to  choose  from  a  list  of  options,  in  which  case  invalid 
answers  will  be  rejected  and  the  question  repeated;  still  others  require 
a  numerical  value  such  as  SNR,  and  any  number  will  be  accepted.  The 
answers  to  some  questions  will  control  the  sequence  of  succeeding 
questions.  All  questions  have  a  default  answer,  so  that  any  question 
which  is  not  given  an  answer  by  the  user  will  nonetheless  be  assigned 
the  default  answer  (usually  a  typical  answer  or  an  answer  which  allow 
the  user  maximun  flexibility  in  succeeding  questions) • 

The  question-answer  section  of  the  program  is  composed  of  three 
parts,  which  are  discussed  in  the  following. 


( t )  The  first  part  requests  some  general  information  about  the 
analysis,  i.e.,  the  type  of  spread-spectrun  technique  to  be  used,  the 
type  of  jamming  present,  and  the  type  of  error-correcting  coding  to  be 
used  ( if  any) . 

( 2 )  The  second  part  requests  the  format  of  the  graphical  output 
desired  by  the  user.  The  user  may  select  both  the  x-axis  variable  and 
the  family  variable.  The  x-axis  variable  is  the  signal  or  jamming 
parameter  that  the  user  wishes  to  vary  in  order  to  evaluate  its  effect 
on  the  bit  or  word  error  rate  of  the  communication  system.  This 
variable  appears  as  the  x-axis  variable  on  the  graphical  plot  output 
produced  by  the  analysis .  Hie  family  variable  is  the  signal  or  jamming 
parameter  which  is  allowed  to  assume  a  number  of  discrete  values 
(maximum  of  9  values)  in  order  to  produce  a  family  of  curves  on  a  single 
plot,  with  each  curve  representing  one  of  the  specified  family  variable 
values.  Hie  user  is  prompted  by  the  program  for  the  x-axis  variable  and 
its  range  of  values,  and  for  the  family  variable  and  its  discrete 
values.  Hie  parameters  that  follow  may  be  specified  as  the  x-axis  or 
family  variables. 

(a)  Number  of  repetitions  of  an  information  bit  when  repetition 
coding  is  used. 

(b)  Number  of  hopping  channels  that  are  receiving  an 
interfering  (i.e.,  jamming)  signal. 

(c)  Signal-bo-noise  ratio  (SNR)  in  decibels. 

(d)  Jamming- to-signal  ratio  (JSR)  in  decibels. 
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(e)  Portion  of  total  spread- spec trim  bandwidth  over  which  an 
interfering  signal  is  present  (i.e  ,  portion  of  band  that  is  jammed) . 
(This  parameter  is  directly  related  to  parameter  (b),  and  as  such, 
parameters  (b)  and  (e)  are  mutually  exclusive.) 

These  five  parameters  are  the  only  ones  that  may  be  varied 
within  a  particular  analysis,  and  only  two  of  these  (the  parameters 
selected  to  be  the  x-axis  and  family  variables)  may  be  varied  at  any  one 
time.  All  other  parameters  (such  as  bandwidths)  are  held  constant 
during  an  analysis. 

(3)  The  third  portion  of  the  question-answer  section  of  the  program 
allows  the  user  to  set  the  values  of  the  constant  parameters  and  the 
parameters  not  selected  as  the  x-axis  or  family  variable.  All  of  these 
remaining  parameters  have  default  values  which  are  used  if  the  user 
declines  to  specify  a  value. 

After  the  question-answer  section  of  the  program  has  been 
completed,  the  user  is  provided  with  a  tabular  synopsis  of  the  parameter 
values  that  will  be  used  in  the  analysis.  The  program  then  runs  the 
analysis  and  produces  the  x-y  plot  containing  the  results  of  the 
analysis;  finally,  the  program  returns  to  the  beginning  of  the  question- 
answer  section  to  allow  the  definition  of  the  parameters  for  a  new 
analysis,  if  desired  by  the  user. 

A  listing  of  the  program  source  code  appears  in  appendix  A. 

Appendix  B  gives  an  example  of  the  interactive  dialogue  between 
the  user  and  the  program  for  a  typical  analysis,  along  with  the  results 
of  the  analysis. 

All  program  variables  of  interest  are  described  in  detail  in 
appendix  C,  and  appendix  D  describes  all  library  subroutines  and 
functions  required  by  the  program.  Finally,  appendix  E  gives  default 
values  for  certain  program  variables.  These  values  are  used  unless  the 
user  specifies  other  values  when  asked. 

The  remainder  of  section  3  is  devoted  to  descriptions  of  the 
internal  structures  of  the  various  routines  used  in  the  program.  All 
the  routines  are  written  in  FORTRAN  IV  and  may  be  executed,  with  minor 
modifications,  on  most  machines  with  sufficient  memory  and  a  FORTRAN  IV 
compiler. 


3. 1  Main  Program 

The  main  program  consists  almost  entirely  of  calls  to 
subroutines  which  perform  the  actual  work  required  for  an  error 
analysis.  The  main  program  first  calls  subroutines  INITP  (which 
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performs  various  initialization  functions)  and  ELT1  (which  zeros  the 
central  processor  unit  timer).  The  primary  program  loop,  starting  at 
ISN  IS,  is  then  entered.  Each  iteration  of  this  loop  constitutes  a 
complete  analysis  of  a  system,  including  plots  and  tabular  output,  The 
loop  begins  by  reinitializing  the  graphics  screen  (ISM  15-17)  then 
displays  the  amount  of  CPU  time  used  during  the  previous  loop  iteration 
(ISN  18-19)  and  reinitializes  the  CPU  timer  (ISN  20).  Then  subroutine 
MODE  is  called  in  ISN  21  to  determine  the  type  of  analysis  to  be 
performed  during  the  current  loop  iteration  (bit  or  word  error  rate, 
type  of  jamming,  type  of  data  modulation,  etc) .  Subroutine  DFAM,  called 
in  ISN  22,  determines  which  parameter  (if  any)  is  to  be  varied  to 
produce  a  family  of  curves  and  the  value  of  the  parameter  to  be  used  for 
each  curve.  Subroutine  DXAX,  called  in  ISN  23,  similarly  determines 
which  parameter  is  to  be  the  x-axis  variable  and  the  range  of  values  for 
the  x-axis.  Subroutine  PARSET  (ISN  24)  then  determines  the  values  of 
all  other  parameters.  Subroutine  TABOUT  (ISN  25)  produces  the  tabular 
output,  including  the  important  characteristics  of  the  analysis  and  the 
values  of  all  input  parameters.  Subroutine  PCALC  then  calculates  the 
bit/word  error  rate  curves  (for  frequency  hopping  and  optionally  for  PN 
modulation,  respectively,  in  ISN  26  and  31)  .  These  curves  are  then 

plotted  on  the  graphics  screen  by  subroutine  FHPLOT  (ISN  27  and  35). 
The  loop  is  then  terminated  in  ISN  37  by  a  transfer  bo  statement  10. 

3. 2  Subroutine  FHPLOT 


Subroutine  FHPLOT  actually  plots  the  analysis  results  on  the 
screen  of  the  graphics  terminal.  The  routine  begins  by  clearing  the 
screen  and  setting  graphics  display  options  concerning  the  actual  plot 
appearance  (ISN  7-17).  The  data  to  be  plotted  are  then  scaled  (ISN  18- 
37)  and  plotted  (ISN  38-56),  and  axis  labels  are  added  (ISN  57-61). 

3.3  Entry  INITP 

This  entry  point  in  subroutine  FHPLOT  initializes  the  plot 
library  when  the  program  is  first  begun.  Axis  labels  defined  in  block 
data  are  processed  for  use  by  the  plot  library  (ISN  66-72).  The  actual 
plotting  area  on  the  terminal  screen  is  then  determined  (ISN  76-86). 
Finally,  the  program  introduction  is  displayed  on  the  screen  if  the  user 
so  desires  (ISN  87-96). 

3.4  Subroutine  MODE 


This  subroutine  determines  the  basic  signal  and  jamming 
characteristics  to  be  used  in  the  analysis.  The  type  of  jamming  to  be 
analyzed  and  some  associated  parameters  are  requested  from  the  user  (ISN 
25-40).  The  coding  parameters  are  then  requested  (ISN  42-65).  Finally, 
some  information  is  requested  concerning  the  type  of  frequency  hopping 
and  type  of  data  modulation  to  be  analyzed  (ISN  69-83). 
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3. 5  Subroutine  DC  AM 


This  routine  determines  the  parameter  to  be  used  as  the  family 
variable.  The  faaiily  variable  is  first  requested  ( ISN  16-31);  then  the 
number  of  family  curves  and  their  corresponding  parameter  values  are 
obtained  from  the  user  ( ISN  33-44). 

3.6  Entry  FAMC 

This  entry  point  in  subroutine  DFAM  automatically  determines 
the  family  parameter  values  when  repetition  coding  is  selected  in 
subroutine  PARSET. 

3.7  Subroutine  DXAX 


This  routine  determines  the  parameter  to  be  used  as  the  x-axis 
variable.  The  x-axis  variable  is  requested  in  ISN  11-22,  then  the  range 
of  values  for  the  x-axis  is  obtained  (ISN  24-30). 

3.8  Entry  XAXC 

This  entry  point  in  subroutine  DXAX  automatically  determines 
the  x-axis  parameter  range  when  repetition  coding  is  selected  in 
subroutine  PARSET. 

3.9  Subroutine  PARSET 


This  routine  obtains  from  the  user  the  values  of  all  parameters 
not  previously  set.  Default  values  are  assigned  in  ISN  17-27.  Defaults 
are  overridden  by  the  user,  if  desired,  in  ISN  28-83.  Hie  use  of 
repetition  coding  is  then  determined  (ISN  84-90)  as  follows:  (1)  if  no 
other  coding  is  specified  in  subroutine  MODE,  repetition  coding  is 
assumed,  and  (2)  the  number  of  repetitions  is  determined  by  the  chip 
bandwidth  (PCM)  divided  by  the  information  bandwidth  (FB).  This  latter 
quotient  is,  by  default,  one,  unless  otherwise  specified  by  FCM  and  FB, 
because  FCM  and  FB  are  equal  by  default.  Thus  simple  bit  error  rate 
(one  repetition)  is  calculated  as  the  default.  Finally,  if  the  nunber 
of  repetitions  is  to  be  the  family  or  x-axis  variable,  then  the 
appropriate  routine  is  called  (ISN  92-94). 

3.10  Subroutine  TABOUT 


This  routine  provides  a  tabular  summary  of  the  important 
parameters  to  be  used  in  the  analysis.  Constants  are  printed  in  ISN  13- 
34.  Family  and  x-axis  parameters  and  their  values  are  then  printed  (ISN 
35-41).  Finally,  values  that  are  in  units  of  decibels,  kilohertz,  or 
megahertz  are  converted  to  their  absolute  values  for  computational  use 
(ISN  42-46) . 
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3.11  Subroutine  PCAIiC 

This  routine  directs  the  calculation  of  the  data  points  on  each 
of  the  family  curves.  The  correct  family  and  x-axis  values  are  assigned 
to  the  family  and  x-axis  variables  for  each  data  point  to  be  computed. 
Then  each  data  point  is  actually  calculated  by  calling  the  appropriate 
function  (FHOP  for  frequency  hopping,  PRPN  for  pseudonoise)  as 
determined  by  the  argument  SUBR,  which  is  assigned  a  value  of  either 
FHOP  or  PRPN  by  the  main  program  based  on  the  type  of  plot  desired  by 
the  user.  The  data  point  is  placed  in  the  array  Y  for  subsequent 
plotting. 

3.12  Function  FHOP 


This  routine  actually  calculates  the  bit  or  word  error  rates 
for  fast  and  slow  frequency  hopping.  Signal  and  jamming  parameters  are 
set  and  checked  in  ISN  10-31.  Intermediate  probabilities  for  use  in  the 
error  rate  evaluation  are  then  calculated  (with  different  probabilities 
for  different  types  of  signals  and  jamming)  in  ISN  33-65.  The  error 
rate  is  then  calculated  for  FSK  modulated  fast-frequency  hopping  with 
narrowband  jamming  in  ISN  66-97,  which  evaluates  equation  (14)  of 
Torrieri.1  Actually,  this  equation  has  seven  nested  sums,  four  of  which 
are  calculated  here,  while  the  inner  three  are  calculated  in  function 
PXKI  (ISN  90).  Error  rates  for  repeater  jamming,  slow  frequency 
hopping,  white  noise  repeater  jamming,  and  coherent  PSK  modulation  are 
determined  in  ISN  98-137. 

3.13  Subroutine  QINT 

This  routine  is  used  by  the  integration  routine  NL9  (see 
function  FHOP,  ISN  52)  to  integrate  a  variation  of  Marcum's  Q  function 
for  use  in  determining  one  of  the  intermediate  probabilities  needed  to 
calculate  error  rates. 

3.14  Function  PXKI 


This  function  is  called  by  function  FHOP  (see  ISN  90  of  FHOP) 
to  evaluate  the  inner  three  sums  of  equation  (14)  of  Torrieri.* 


*/?.  J.  Torrieri,  Frequency  Hopping  in  a  Jamming  Environment ,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCH-78-2 
( December  1978 )  . 
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3.15  Function  PRPN 


This  function  la  called  by  subroutine  PCAIC  to  calculate  the 
error  rates  for  pseudonoise  modulated  signals.  This  routine  evaluates 
equations  (54)  and  (55)  of  Tbrrieri.2 

3. 16  Subroutine  FREAD 


This  routine  reads  one  real  number  from  the  terminal  without 
performing  a  line  feed  before  the  read  takes  place.  This  routine  may  be 
.replaced  by  a  READ  statement,  but  this  will  not  allow  the  user's  answer 
to  a  question  to  follow  the  question  on  the  same  line  on  the  terminal 
screen . 

3.17  Subroutine  I READ 


This  routine  reads  one  integer  number  from  the  terminal  without 
performing  a  line  feed  before  the  read  takes  place.  This  routine  may  be 
replaced  by  a  READ  statement • 

3.18  Function  LOGCOM 


This  routine  computes  the  natural  logarithm  of  ( I ,  J)  where  I 
and  J  are  the  arguments  of  the  function  and 


(t,J) 


II 

J!  (I  -  J)! 


3.19  Function  FACT 

This  routine  calculates  I!  ,  where  I  is  the  argument  of  the 

function. 

3.20  Block  Data 

This  routine  initializes  the  x-  and  y-axis  label  alphanumeric 
values  for  use  by  the  plotting  routines. 


2D.  J.  Torrieri,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  Report  CH/CCM-79-9 
( December  1979)  . 
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4.  SUMMARY 
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A  computer  program  has  been  developed  which  analyzes  and  compares 
the  performance  of  two  types  of  spread-spectrum  techniques  in 
interfering  or  jamming  environments.  The  program  does  not  consider 
synchronization  problems  or  signal  and  jamming  propagation  effects. 
Instead,  the  program  assumes  that  the  transmitter  and  receiver  are  in 
synchronization  and  that  there  is  a  defined  signal-to-noise  ratio  and 
jamming-to-signal  ratio  at  the  receiver.  The  program  can  give 
considerable  insight  into  the  performance  increases  and  relative  values 
that  might  be  achieved  by  the  incorporation  of  spread-spectrum 
techniques  into  existing  or  postulated  communication  systems. 
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APPENDIX  A. — PROGRAM  LISTING 
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APPENDIX  A — PROGRAM  LISTING 

This  appendix  contains  a  listing  of  the  spread-spectrum  analysis 
program.  Further  details  concerning  the  use  of  the  program  as  well  as 
descriptions  of  subroutines  and  variables  used  in  the  program  are  given 
in  appendices  B,  C,  D  and  E. 
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APPENDIX  B.  —SAMPLE  PROGRAM  DIALOGUE 

The  program  begins  by  asking  if  the  user  wishes  an  overview  of  the 
program  (see  fig.  B-1 ) .  If  the  user  answers  "y,"  then  the  brief  program 
description  is  displayed  before  the  program  proceeds  bo  the  question- 
answer  part  of  the  program.  If  the  user  answers  "n,“  then  the 
description  is  not  displayed.  (All  user  responses  in  the  dialogue  are 
underlined . ) 


TYPE  Y  FOR  OVERVIEW  OF  PROGRAM, 

THIS  PROGRAM  ESTIMATES  THE  BIT  ERROR  RATE  FOR  FSK  FREQUENCY 
HOPPING  AND  PN  COMMUNICATION  SYSTEMS.  NOISE  AND  JAMMING  ARE 
CONSIDERED  IN  THE  MODEL.  THIS  PROGRAM  MUST  BE  RUN  FROM  A  TEKTRONIX 
GRAPHICS  TERMINAL  OR  EQUIVALENT . 

TO  HALT  EXECUTION  OF  THIS  PROGRAM,  TYPE  t  IN  RESPONSE 
TO  ANY  REQUEST  FOR  NUMERICAL  INPUT. 

TYPE  RETURN  TO  CONTINUE  _ 

Figure  B-1.  Sample  brief  program  description. 


The  question-answer  program  begins  by  displaying  the  amount  of 
Central  Processing  Unit  (CPU)  time  used  by  the  previous  analysis  (if 
there  was  one)  ,  and  then  begins  asking  questions  in  order  to  define  the 
analysis  parameters  (see  fig.  B-2).  The  type  of  jamming  to  be  used  is 
first  determined.  Next,  the  user  is  asked,  "Constant  jamming  power?" 
If  the  answer  to  this  question  is  yes,  then  the  total  amount  of  jamming 
power  is  held  constant  in  any  analyses  which  vary  the  number  of  jammed 
channels  or  the  portion  of  the  band  to  be  jammed  (either  as  the  x-axi3 
or  the  family  variable).  Thus,  when  varying  the  number  of  jammed 
channels  or  the  portion  of  the  band  to  be  jammed,  the  program  dilutes  or 
concentrates  the  jamming  power  density  accordingly,  keeping  the  total 
jamming-to-signal  ration  ( JSR)  a  constant  over  the  total  bandwidth.  If 
the  user  answers  "no"  to  this  question,  the  jamming  power  density  is 
held  constant,  and  an  increase  or  decrease  in  the  number  of  jammed 
channels  or  portion  of  the  band  jammed  will  cause  a  corresponding 
increase  or  decrease  in  the  total  JSR  over  the  total  bandwidth. 

The  user  is  then  asked  to  specify  the  data-coding  parameters  to  be 
used.  Requests  are  made  for  word  length  (the  number  of  bits  per  word 
before  coding),  chips  per  word  (where  a  "chip"  is  a  coded  bit  and,  thus, 
chips  per  word  is  the  number  of  bits  per  word  after  coding)  ,  and  chip 
error  threshold  (the  minimum  number  of  coded  bits  which  must  be  in  error 
in  order  to  prevent  the  coding  from  correcting  these  errors) .  If  the 
word  length  is  specified  to  be  one,  then  bit  error  rates  are  calculated 
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and  the  remaining  questions  concerning  coding  are  skipped.  (It  should 
be  noted  that  bit  error  rates  may  be  calculated  for  either  uncoded  data 
or  repetition-coded  data.  Repetition  coding  and  the  method  for 
specifying  its  use  in  the  program  are  explained  in  the  description  of 
the  third  part  of  the  question-answer  portion,  below).  *nie  user  is  then 
asked  if  frequency  hopping  and  peeudonoise  curves  are  to  be  plotted  on 
the  sa.T*e  set  of  axes  or  on  two  different  sets.  Following  this,  the  user 
is  asked  "Do  you  want  expanded  integrals?"  A  "no"  answer  causes  certain 
integrals  used  in  the  calculation  of  the  bit  error  rates  to  be  estimated 
rather  than  calculated.  This  estimation  is  used  primarily  in  special 
cases  either  where  the  estimation  is  known  to  be  valid  or  where  the 
evaluation  of  the  integrals  is  too  costly  in  CPU  time.  Next,  the  user 
is  asked  if  slow  frequency  hopping  is  to  be  analyzed  instead  of  fast 
hopping.  If  slow  hopping  is  selected,  the  user  is  finally  asked  if 
coherent  phase  shift  keying  (PSK)  data  modulation  is  to  be  analyzed 
instead  of  frequency  shift  keying  (FSK)  modulation  (the  default). 


TINE  =  0:00.01'’ 

JANNING  OPTIONS: 

1  NARROU  BAND 

2  PARTIAL  BAND 

3  REPEATER 
JAPIMING  TYPE:  _i_ 

CONSTANT  JANNING  POUER  (V  OR  N)7_y_ 

IF  UORD  LENGTH  IS  ONE  BIT,  BIT  ERROR  RATE  IS  CALCULATED. 
OTHERUISE  UORD  ERROR  RATE  IS  CALCULATED. 

UORD  LENGTH  IN  BITS:  A 
CHIPS  PER  UORD:  4 
CHIP  ERROR  THRESHOLD: 

HOPPING,  PN  CURVES  ON  SANE  PLOT?^ 

DO  YOU  UANT  EXPANDED  INTEGRALS  (Y  OR  N )?  y 

SLOU  FREQ  HOPPING7  n 


Figure  B-2.  Question-answer  portion — part  1. 
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Ttie  second  part  of  the  question-answer  portion  of  the  program  is 
shown  in  figure  B-3.  The  user  is  asked  to  specify  the  parameters  to  be 
used  for  the  x-axis  and  family,  and  their  values  or  ranges.  The  x-axis 
increment  is  the  interval  between  adjacent  points  evaluated  along  the  x- 
axis  in  producing  the  data  curves.  If  the  x-axis  limits  and  increments 
indicate  that  more  than  100  points  are  to  be  evaluated,  then  the  x-axis 
upper  limit  is  changed  so  that  only  100  points  are  evaluated. 


VALID  FAHILY  AND  X-AXIS  VARIABLES  ARE: 

6  NONE  (FAMILY  ONLY) 

1  NBR  OF  CHIPS  PER  BIT 

2  NBR  OF  JAMMED  CHANNELS 

3  SIGNAI-T0-N0ISE  RATIO  ( FSK ) 

A  JAMMING-TO-SIGNAL  RATIO 


FAMILY  (t  TO  STOP)?  A 

NBR  OF  FAMILY  CURVES:  _JL 
JSR  IN  DB? 

JSR  IN  DB?  _2fi_ 

JSR  IN  DB?  2$L 


X-AXIS?  JL 

X-AXIS  START?  A3_ 

X-AXIS  END? 

X-AXIS  INCR 

Figure  B-3.  Question-answer  portion — part  2. 

The  third  part  of  the  question-answer  portion  is  shown  in  figure  B- 
4.  All  variables  not  previously  assigned  values  are  displayed,  along 
with  their  default  values.  If  all  these  default  values  are 

satisfactory,  then  the  user  need  only  respond  "yes"  to  the  question, 
"Are  defaults  OK?”  If  the  user  answers  "no,”  then  he  will  be  prompted 
for  values  for  each  of  the  parameters  that  may  legally  be  specified  at 
this  time.  If  the  user  types  only  a  return  in  response  to  any  prompt, 
the  default  value  for  that  parameter  will  be  assumed. 

The  parameter  values  specified  or  defaulted  in  the  third  question- 
answer  portion  are  used  to  determine  certain  other  parameter  values. 
The  n.imher  of  hopping  channels  available  to  the  frequency  hopping  system 
to  be  analyzed  is  letermined  by 


_2fi_ 


NTH 


A 


(  FW 


NBPW)  /(  NBPC 


FB) 
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***  VARIABLE  DEFAULTS  ttt 

NBR  OF  JAMMED  CHANNELS  ( NJAM )  *  1 

PORTION  JAMMED  (UPBJ)  •  0.2500 
INFO  BIT  BANDWIDTH  (FB)  *  25.0  KHZ 

PN  PROCESSING  GAIN  <G)  »  30.0  DB 

HOPPING  BANDWIDTH  (BU)  «  25.0  MHZ 

CHIP  BANDUIDTH  (FCM )  -  25.0  KHZ 

HOPPING  OVERLAP  FACTOR  (A)  «  1.0000 
COMM-TO-JAM  DIFF  FREQ  FACTOR  (B)  -  1.0000 
PN  RCVR  LOSS  FACTOR  (D)  •  1.0000 

ARE  DEFAULTS  OK  (V  OR  N)?_n_ 

NBR  OF  JAMMED  CHANNELS: — 

INFO  BIT  BANDUIDTH  (FB  IN  KHZ): — 

PN  PROCESSING  GAIN  (G  IN  DB): — 

HOPPING  BANDUIDTH  (BU  IN  MHZ): — 

HOPPING  OVERLAP  FACTOR  (A):— 

COMM-TO-JAM  DIFF  FREQ  FACTOR  (B): — 

PN  RCVR  LOSS  FACTOR  (D)  :  — 


Figure  B-4.  Question-answer  portion — part  3. 


where  NBPW  and  NBPC  are  the  number  of  information  bits  per  word  and  the 
number  of  coded  bits  per  word,  respectively,  and  where  A  (hopping 
channel  overlap  factor),  BW  (hopping  bandwidth),  and  FB  (information  bit 
bandwidth)  are  set  or  defaulted  in  the  third  question-answer  portion. 
In  addition,  repetition  coding  may  also  be  specified  in  this  portion. 
Repetition  coding  is  used  to  improve  bit  error  rates  by  transmitting 
each  information  bit  a  number  of  times  with  a  change  in  frequency  (  in 
hopping  systems)  between  each  repetition,  and  using  majority  logic  in 
the  receiver  to  determine  the  true  value  of  the  information  bit.  The 
use  of  repetition  coding  is  controlled  by  the  values  of  FCM  (coded  bit 
bandwidth)  and  FB  ( information  bit  bandwidth)  when  bit  error  rates 
(rather  than  word  error  rates)  have  been  previously  specified.  If  PCM 
and  FB  are  equal  (which  is  the  default)  then  no  repetition  coding  is 
used,  and  simple  bit  error  rates  are  calculated.  If  PCM  and  FB  are  not 
equal,  then  repetition  coding  is  assumed,  with  the  number  of  repetitions 
specified  by 


NBPC  =  FCM/  FB . 
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NBPC  is  an  integer  and  is  therefore  truncated  to  an  integar  value  if  the 
ratio  of  PCM  and  FB  is  not  a  whole  number.  If  FCM  is  smaller  than  FB, 
then  NBPC  is  set  to  one  and  the  simple  bit  error  rate  is  calculated. 

Figure  B-5  shows  the  tabular  output  produced  by  the  program  and 
summarizes  the  important  analysis  parameters  and  their  values.  The 
tabular  output  will  remain  on  the  display  screen  until  the  user  hits  the 
return  key  and  allows  the  program  to  continue.  Figure  B-6  shows  the 
analysis  results  in  the  form  of  a  plot  of  the  requested  data.  The  solid 
line  on  the  plot  is  actually  an  overlapping  of  the  three  family  curves 
computed  for  frequency  hopping.  This  overlapping  occurred  because  the 
variation  of  the  family  variable  { jamming-to-signal  ratio)  had  no  effect 
on  the  performance  of  the  frequency  hopping  signal .  The  three  family 
curves  for  pseudonoise  modulation  appear  as  dashed  lines  in  order  to 
differentiate  them  from  the  frequency  hopping  curves  when  both  sets  of 
curves  are  plotted  on  the  same  set  of  axes.  The  curves  are  not 
automatically  labeled  by  the  program  with  their  corresponding  family 
variable  values  because  of  tbe  difficulty  in  determining  where  such 
labeling  should  be  placed.  Thus,  the  user  must  watch  the  display  screen 
as  the  plot  is  generated  in  order  to  determine  the  correct  family 
variable  value  for  each  curve  (the  curves  are  drawn  in  the  order  in 
which  the  family  values  were  specified  by  the  user  in  the  question- 
answer  portion  of  figure  B-3 ) .  Hie  plot  will  remain  on  the  display 
screen  until  the  user  hits  the  return  key.  The  program  will  then  return 
to  the  first  part  of  the  question-answer  portion  (fig.  B-2 )  to  allow  the 
user  to  begin  a  new  analysis. 


HOPPING  RATE  =  FAST 

JAMMER  TYPE  =  NARROW  BAND  ERROR  RATE  *  UORD 

4  BITS  PER  UORD  1  URONG  FOR  UORD  ERROR 

NBR  OF  CHIPS  PER  UORD  *  4 

NBR  OF  JAMMED  CHANNELS  =  1 

INFO  BIT  BANDUIDTH  *  25.0  KHZ 

PN  PROCESSING  GAIN  «  30.0  DB 

HOPPING  BANDUIDTH  *  25.0  MHZ 

HOPPING  OUERLAP  FACTOR  =  1.00 
COMM-TO-JAM  DIFF  FREQ  FACTOR  *  1.00 
PN  RCUR  LOSS  FACTOR  «  1.00 

FAMILY:  JAMMING-TO-SIGNAL  RATIO  (DB) 

UALUES  UERE:  10.0,  20.0,  30.0, 

X-AXIS:  SIGNAL/NOISE  PER  UORD  (DB) 

START:  13.0  END:  26.0 
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APPENDIX  C. --COMMON  VARIABLES 

All  program  variables  of  interest  to  the  user  are  described  in 
detail  in  the  following  appendix. 

COMMON  BLOCK  PLOT 

X(100)  Real  array  containing  up  to  100  x  coordinates  (i.e., 

x-axis  values)  to  be  plotted.  The  exact  number  of 
points  in  X  is  given  by  the  variable  NPTSX  in  COMMON 
Block  XAXPAR.  The  values  in  X  are  set  in  subroutine 
DXAX  (ISN  40-42)  and  are  used  in  subroutine  PCALC  to 
calculate  the  error  probability  curves. 

Y(  100,9)  Real  array  containing  up  to  nine  family  curves  with 

up  to  100  y  coordinates  (i.e.,  y-axis  values)  to  be 
plotted .  The  exact  number  of  family  curves  to  be 
plotted  is  given  by  the  variable  NFAM  in  COMMON  block 
FAMPAR.  The  exact  number  of  points  in  each  curve  is 
given  by  the  variable  NPTSX  in  COMMON  block  XAXPAR. 
Hie  values  in  Y  are  set  in  subroutine  PCALC.  This 
variable  corresponds  to  the  probability  of  error,  P  , 
in  equation  (14)  of  Torrieri.1 


COMMON  BLOCK  INPAR 

NJAM  Integer  variable  whose  value  is  the  number  of 

distinct  frequency  hopping  channels  which  are  being 
jammed  by  CW  (tone)  or  wide-band  noise  jammers.  This 
variable  may  be  the  x-axis  variable,  the  family 
variable,  or  a  constant.  NJAM  must  be  less  than  the 
total  number  of  hopping  channels,  and  must  not  be 
less  than  zero.  This  variable  corresponds  to  the 
number  of  channels  jammed,  j,  in  Torrieri.' 

BETA  1  Real  variable  whose  value  is  the  received  signal-to- 

noise  ratio  (SNR)  per  coded  word  in  decibels.  It  is 
converted  to  SNR  per  coded  bit  in  subroutine  TABOUT 
(ISN  42)  or  subroutine  PCALC  (ISN  15  or  26),  and  is 
converted  to  SNR  per  uncoded  (i.e.,  information)  bit 
in  subroutine  FHOP  (ISN  26).  This  variable  may  be 
the  x-axis  variable,  the  family  variable,  or  a 


1D.  J.  Torrieri,  Frequency  Hopping  in  a  Jamming  Environment,  V.S> 
Army  Development  and  Read iness  Command,  DARCOM  Report  CH/CCM-78-2 
(December  1978). 
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constant.  BETA1  corresponds  to  wRg/Ntu  for  frequency 
hopping  (Torrieri,1  p  21)  and  to  R /N0Bm  for  W 
modulation.  (Rs/N0Bm  is  related  to  N0/Eb  in  equation 
(55)  of  Torrieri2  by  V*o  =  (VNoBm)BmTm'  BmTm 

corresponds  to  program  variable  D] 

GAMMA  Real  variable  whose  value  is  the  received  jamming-to- 

signal  ratio  (JSR)  in  decibels.  For  frequency 
hopping ,  GAMMA  may  be  either  the  JSR  per  frequency 
hopping  channel  or  the  JSR  over  a  portion  or  all  of 
the  entire  hopping  bandwidth  (i.e.,  the  jammer  power 
is  spread  out  over  a  number  of  hopping  channels) , 
determined  by  the  value  of  the  variable  ICON  in  this 
COMMON  block.  The  value  .of  ICON  is  set  by  the  use 
in  subroutine  MODE.  For  modulation,  GAMMA  is  the 
total  JSR  due  to  all  jamming  signals  in  the  PN 
bandwidth.  GAMMA  may  be  the  x-axis  variable,  the 
family  variable,  or  a  constant.  GAMMA  corresponds  to 
Rj/Rs  in  Torrieri.1,2 

FB  Real  variable  whose  value  is  the  information 

(uncoded)  bit  bandwidth  in  kilohertz  (i.e.,  the 
information  rate  in  kilobits  per  second  of  the 
transmitted  bit  stream  if  it  were  uncoded)  .  Ibis 
variable  is  a  constant  set  by  the  user  in  subroutine 
PARSET,  and  corresponds  to  2ffa  in  equation  (3)  and  ^ 
in  equation  (48)  of  Torrieri.1 

BW  Real  variable  whose  value  is  the  total  bandwidth  in 

megahertz  available  for  frequency  hopping 

communications.  Ibis  variable  is  a  constant  set  by 
the  user  in  subroutine  PARSET,  and  corresponds  to  Bo 
in  equation  (3)  of  Tbrrieri.1 

FCM  Real  variable  whose  value  is  the  coded  bit  rate  of 

the  transmitted  bit  stream  (i.e.,  the  coded  bit 
bandwidth)  in  kilohertz  for  repetition  coding  only. 
The  number  of  coded  bits  per  uncoded  bit  when 
repetition  coding  is  used  is  calculated  as  NBPC  = 
FCM/FB,  where  NBPC  and  FB  are  as  specified  in  this 
COMMON  block  description.  This  variable  ( FCM)  is  a 
constant  set  by  the  user  in  subroutine  PARSET,  and 
corresponds  to  2fc  in  equation  (3)  of  Torrieri.1 

1 D .  J.  Torrieri,  Frequency  Hopping  in  a  Jamming  Environment,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
(December  1979). 

2D.  J.  Torrieri,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Read iness  Command,  Report  CM/CCM-79-9 
(December  1979). 
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A  Real  variable  whose  value  is  a  parameter  which 

accounts  for  any  hopping  channel  overlap  or 
separation  (A  =  1  for  adjacent  channels,  A  >  1  for 
overlapping  channels,  A  <  1  for  separated 

channels).  This  variable  is  a  constant  set  by  the 
user  in  subroutine  PARSET,  and  corresponds  to  "a"  in 
equation  (3)  of  Tbrrieri.* 

8  Real  variable  whose  value  is  a  parameter  expressing 

the  type  of  jamming  against  TO  systems.  Hie  worst 
type  of  jamming  (with  respect  to  the  TO  system)  is 
continuous  wave  jamming  at  the  center  frequency  of 
the  TO  transmission;  in  this  case,  B  =  2.  For  other 
types  of  jamming,  B  <  2.  This  variable  is  a  constant 
whose  value  is  set  by  the  user  in  subroutine  PARSET, 
and  corresponds  to  "b"  in  equation  (55)  of  Torrieri.^ 

D  Real  variable  whose  value  is  a  parameter  expressing 

the  receiver  losses  in  a  TO  system.  The  value  of  D 
typically  varies  between  1  and  2,  with  2  being  the 
worst  case  (with  respect  to  the  TO  receiver)  .  this 
variable  is  a  constant  whose  value  is  set  by  the  user 
in  subroutine  PARSET,  and  corresponds  to  "W 
( time- bandwidth  product)  in  the  expression  for  N^E^ 
as  described  in  this  section  as  part  of  the 
discussion  of  variable  BETA1. 

NBPC  Integer  variable  whose  value  is  the  number  of  "chips” 

(coded  bits)  per  word.  For  example,  if  (7,4)  coding 
is  used  (where  4  information  bits  are  coded  into  7 
transmitted  bits) ,  NBPC  =  7.  If  coding  is  not 

employed,  NBPC  =  NBPW  (see  the  description  of  NBPW  in 
this  COMMON  block).  If  bit  error  rate  is  desired, 
NBPC  =  NBPW  =1.  If  repetition  coding  is  used,  NBPC 
is  the  number  of  times  each  bit  is  repeated.  This 
variable  is  a  constant  (except  in  repetition  coding, 
where  it  may  be  the  family  or  x-axis  variable)  set  by 
the  user  in  subroutine  MODE,  and  corresponds  to  "C" 
in  equations  (1)  and  (50)  of  Torrieri.1 2 


1D.  J.  Torrieri ,  Frequency  Hopping  in  a  Jamming  Environment ,  U.S . 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCH-7 8-2 
(December  1978). 

2D.  J.  Torrieri ,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  Report  CH/CCM-79-9 
(December  1979). 
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IWN  Integer  variable  whose  value  specifies  the  type  of 

jammer  modulation  used  by  a  repeater  jammer.  If  IWN 
=  1,  then  white-noise  modulation  is  used.  If  IWN  = 
0,  then  no  modulation  is  used  (i.e.,  the  jammer 
merely  retransmits  its  received  signal) .  This 
variable  is  a  flag  set  by  the  user  in  subroutine 

MODE. 

IFP  Integer  variable  whose  value  specifies  the  use  of 

coherent  phase- shift-keyed  (PSK)  data  modulation 

rather  than  frequency-shift-keyed  (FSK)  modulation 

for  frequency  hopping.  If  IFP  =  1,  then  coherent  PSK 
modulation  is  used;  if  IFP  =  0,  FSK  is  used.  This 
variable  is  a  flag  set  by  the  user  in  subroutine 

MODE. 

G  Real  variable  whose  value  is  the  processing  gain  in 

decibels  of  the  PN  system.  This  variable  is  a 
constant  set  by  the  user  in  subroutine  PARSET,  and 
corresponds  to  G  in  equation  (55)  of  Torrieri.2 

NBPW  Integer  variable  whose  value  is  the  number  of 

information  (i.e.,  uncoded)  bits  per  word  in  the  data 
transmission.  For  example,  if  (7,4)  coding  is  used, 
NBPW  =  4.  If  NBHtf  =  1,  then  bit  error  rates  are 

computed  in  the  program;  otherwise,  word  error  rates 
are  computed,  and  the  user  is  issued  iequests  for 
coded  word  length  and  coding  parameters.  This 
variable  is  a  constant  set  by  the  user  in  subroutine 
MODE,  and  corresponds  to  "w"  in  equation  (  1  )  of 
Torrieri . * 

NERR  Integer  variable  whose  value  is  the  smallest  number 

of  coded  bits  per  word  which  must  be  received  in 
error  in  order  to  cause  a  word  error  at  the 
receiver.  For  example,  (7,4)  coding  allows  correct 
determination  of  the  four  information  bits  if  no  more 
than  one  of  the  seven  coded  bits  is  received  in 
error .  Thus ,  at  least  two  bit  errors  must  occur  in  a 
word  in  order  for  a  word  error  to  occur;  NERR  =  2  in 
this  case.  This  variable  is  a  constant  set  by  the 
user  in  subroutine  MODE,  and  corresponds  to  "r"  in 
equation  (2)  and  (50)  of  Torrieri.1 


1 D .  J.  Torrieri ,  Frequency  Hopping  in  a  Jamming  Environment,  U.S. 
Army  Development  and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
(December  1978). 

2D.  J.  Torrieri,  Pseudonoise  Spread-Spectrum  Systems  in  Communication 
Warfare,  U.S.  Army  Development  and  Readiness  Command,  Report  CM/CCM-79-9 
(December  1979). 
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JTP  Integer  variable  whose  value  specifies  the  type  of 

jamming  employed.  The  possible  values  of  JTP  are  JTP 
=  1  (narrowband,  or  continuous  wave,  jamming);  JTP  = 
2  (partial  band  or  wideband  jamming);  JTP  =  3 
(repeater  jamming).  This  variable  is  a  flag  set  by 
the  user  in  subroutine  MODE. 

UPBJ  Real  variable  whose  value  is  the  portion  of  the  total 

hopping  or  PN  bandwidth  ( EW  or  BS  in  this  COMMON 
block)  that  is  being  jammed  by  one  or  more  wideband 
noise  jammers.  This  variable  is  functionally  related 
to  the  variable  NJAM  (see  the  description  of  NJAM  in 
this  COMMON  block).  Thus,  only  one  of  the  two  may  be 
set  by  the  user — the  other  is  determined  from  the 
value  of  the  one  that  is  set  (see  lines  ISN  12  and  14 
in  subroutine  FHOP) .  This  variable  may  be  the  x-axis 
variable,  the  family  variable,  or  a  constant.  The 
value  of  UPBJ  must  not  be  less  them  zero  nor  greater 
than  one.  This  variable  corresponds  to  **u"  in 
equation  (46)  of  Torrieri.* 

IEST  Integer  variable  whose  value  specifies  whether  or  not 

an  approximation  to  equation  (14)  in  Torrieri*  is 

used.  If  IEST  =  1,  the  approximation  is  not  used, 

and  equation  (  14)  is  calculated  in  its  entirety.  If 

IEST  =  0,  then  the  summation  over  q  in  equation  (14) 
is  not  performed  (i.e.,  q  is  set  to  zero),  and  the 
value  of  P3  in  equation  (13)  is  set  to  one.  This 
approximation  considerably  speeds  up  program 
execution  since  S2  in  equation  (13)  need  not  be 

computed  (the  value  of  S2  is  determined  by  equation 
(28)  and  involves  the  time-consuming  integration  of  a 
highly  complex  function).  However,  the  accuracy  of 
the  results  suffers,  particularly  where  large  numbers 
of  jammers  are  present  or  a  significant  portion  of 
the  band  is  being  jammed.  The  approximation  is 
normally  not  used  unless  it  is  certain  that  no 
significant  computation  errors  will  occur.  This 
variable  is  a  flag  set  by  the  user  in  subroutine 
MODE. 


1D.  J.  Torrieri,  Frequency  Hopping  in  a  Jamming  Environment ,  U.S. 
Army  Development  and  Readiness  Command ,  DARCOM  Report  CH/CCM-78-2 
(December  1978). 
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NCH 


ICON 


I  SAME 


ISLOW 


*D.  J.  Tocrieri 
Army  Development 
( December  1978). 


Integer  variable  whose  value  is  the  number  of  hopping 
channels  available  to  the  frequency  hopping  system. 
This  variable  is  calculated  using  equation  (3)  in 
Torrieri*  (see  ISN  11  in  subroutine  FHOP) ,  and 
corresponds  to  the  variable  M  in  that  equation. 

Integer  variable  whose  value  specifies  whether  or  not 
total  jamming  power  is  to  remain  constant  (spread 
over  either  a  portion  of  the  band  or  a  number  of 
continuous  wave  jammers)  .  If  ICON  =  0,  total  jamming 
power '  is  not  constant,  and  as  the  number  of  jammers 
(or  portion  of  the  band  jammed)  increases,  the  amount 
of  jamming  power  increases  proportionally.  If  ICON  = 
1,  then  changes  in  the  number  or  bandwidth  of  the 
jammers  does  not  change  the  total  jamming  power 
(i.e.,  more  jammers  result  in  less  power  per  jammer; 
thus,  jamming  power  is  "diluted").  This  latter  case 
is  used  to  determine  optimum  jamming  strategies  when 
the  jammers  are  limited  in  total  power  available 
(this  is  the  case  in  the  real  world).  The  variable 
ICON  is  a  flag  set  by  the  user  in  subroutine  MODE. 

Integer  variable  whose  value  determines  whether  or 
not  frequency  hopping  and  PN  error  probability  curves 
are  to  be  plotted  on  the  same  plot  for  comparison. 
If  ISAME  =  0,  then  frequency  hopping  curves  are 
plotted  separately  from  the  PN  curves.  If  ISAME  =  1, 
then  both  sets  of  curves  are  plotted  on  the  same 

plot.  This  variable  is  a  flag  set  by  the  user  in 
subroutine  MODE. 

Integer  variable  vdiose  value  specifies  either  fast  or 
slow  frequency  hopping.  If  ISLCW  =  0,  then  the 

equations  for  fast  frequency  hopping  (see  section  2 
of  Torrieri*)  are  used.  If  ISLCW  =  1,  then  the 

equations  for  slow  frequency  hopping  (section  3, 

reference  1)  are  used.  This  variable  is  a  flag  set 
by  the  user  in  subroutine  MODE. 


,  Frequency  Hopping  in  a  Jamming  Environment ,  U.S. 
and  Readiness  Command,  DARCOM  Report  CM/CCM-78-2 
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COMMON 


IF  AM 


NFAM 


FAMV 


COMMON 

IXAX 


BLOCK  FAMPAR 


Integer  variable  whose  value  determines  which  program 
variable  is  to  be  varied  in  order  to  produce  a  family 
(i.e.,  a  series)  of  curves.  The  allowable  family 
variables  and  their  corresponding  values  of  IFAM  are 
none  (IFAM  =  0),  number  of  chips  per  bit  (NBPC,  for 
repetition  coding  only)  (IFAM  =  1),  number  of  jammed 
channels  (NJAM,  cw  jamming  only)  (IFAM  =  2),  signal- 
to-noise  ratio  (BETA1)  per  word  (IFAM  =3),  jamming- 
to-signal  ratio  (GAMMA)  (IFAM  =  4),  and  portion  of 
band  jammed  (UPBJ,  for  partial  band  jamming  only) 
(IFAM  =  5).  13118  variable  is  a  flag  set  by  the  user 
in  subroutine  DFAM. 

Integer  variable  whose  value  is  the  number  of  family 
curves  to  be  calculated.  NFAM  may  be  no  greater  than 
nine.  This  variable  is  a  constant  set  by  the  user  in 
subroutine  DFAM. 

Real  array  containing  up  to  nine  family  variable 
parameter  values.  The  exact  number  of  parameter 
values  in  FAMV  is  given  by  the  variable  NFAM.  A 
separate  error  probability  curve  is  generated  for 
each  of  the  NFAM  values  in  FAMV.  The  parameter  to  be 
assigned  the  values  in  FAMV  iB  given  by  the  variable 
IFAM.  The  actual  assignment  of  the  values  to  the 
parameter  occurs  in  ISN  11-19  of  subroutine  PCALC. 
The  values  of  FAMV  are  set  by  the  user  in  subroutine 
DFAM. 


BLOCK  XAXPAR 


Integer  variable  whose  value  determines  which  program 
variable  is  to  be  the  x-axis  variable.  The  allowable 
x-axis  variables  and  their  corresponding  values  of 
IXAX  are  the  same  as  for  IFAM  in  COMMON  block  FAMPAR 
except  that  IXAX  =  0  (no  x-axis)  is  not  allowed.  Ihe 
variable  selected  to  be  the  x-axis  by  IXAX  is  the 
variable  whose  value  is  varied  to  calculate  each 
point  on  an  error  probability  curve.  The  values  of 
the  x-axis  variable  to  be  used  by  the  program  are 
calculated  using  the  variables  XAXS,  XAXE,  and  XAXI 
in  this  COMMON  block.  The  actual  values  are  then 
placed  in  the  array  X  in  COMMON  block  PLOT;  there  are 
NPTSX  (see  this  COMMON  block)  of  these  values.  TTiis 
variable  is  a  flag  set  by  the  user  in  subroutine 
DXAX. 


65 


APPENDIX  C 


XAXI 


NPTSX 


COMMON  BLOCK 


LAB 


LAB  2 


LAB  3 


LABY 


>*• 


Real  variable  whose  value  is  the  desired  increment 
between  successive  x-axis  variables.  Thus,  the  nth 
x-axis  value  (  in  array  X)  is  calculated  by  adding 
XAXI  to  the  (n-1)t^1  x-axis  value.  It  should  be  noted 
that  since  successive  x-axis  value  are  calculated  by 
adding  the  increment  to  the  previous  x-axis  value, 
the  number  of  x-axis  values  will  exceed  100  (the 
maximum  allowed)  if  XAXI  is  less  than  one  one- 
hundredth  of  the  span  of  the  x-axis  (XAXE-XAXS).  In 
this  case  the  number  of  x-axis  values  is  truncated  to 
100,  and  the  x-axis  will  not  reach  XAXE.  The 
variable  XAXI  is  a  constant  set  by  the  user  in 
subroutine  DXAX. 

Integer  variable  whose  value  is  the  number  of  x-axis 
values  that  are  in  the  array  X  (in  COMMON  block  PLOT) 
to  be  plotted.  This  variable  must  not  be  less  than 
one  nor  greater  than  100.  The  value  of  NPTSX  is 
calculated  in  ISN  35  of  subroutine  DXAX,  using  XAXS, 
XAXE,  and  XAXI. 


LABELS 


Integer  array  dimensioned  to  35  words  and  containing 
the  five  titles  for  the  x-axis  ( see  description  of 
IXAX  in  COMMON  block  XAXPAR) .  Each  title  is  seven 
words  long,  and  the  characters  are  packed  four  to  a 
word.  The  titles  are  placed  in  the  array  at  compile 
time  in  the  BLOCK  DATA  subroutine  (ISN  5).  The  title 
characters  are  initially  in  EBCDIC  but  are  translated 
to  ASCII  in  ISN  69  of  entry  point  INITP  (subroutine 
FHPLOT ) . 

This  array  contains  an  exact  copy  of  the  contents  of 
array  LAB  before  LAB  is  translated  to  ASCII.  LAB2  is 
used  in  subroutine  TABOUT  to  display  the  family  and 
x-axis  variables  to  the  user. 

Integer  array  dimensiond  to  140  and  containing  the 
unpacked  x-axis  titles.  The  titles  are  unpacked  from 
four  to  a  word  (in  array  LAB)  to  one  to  a  word  (LAB 3) 
in  ISN  7  of  entry  point  INITP  (subroutine  FHPLOT). 
Unpacking  is  performed  after  translation  from  EBCDIC 
to  ASCII  character  codes  has  taken  place. 

Integer  array  containing  the  four-word  y-axis 
title.  The  title  is  packed  four  characters  to  a 
word,  and  is  initialy  placed  in  LABY  in  the  BLOCK 
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DATA  subroutine  ( ISN  6).  The  characters  are 
initially  in  ffiCDIC  but  are  translated  to  ASCII  in 
ISN  71  of  entry  point  INITP  (subroutine  FHPLOT) . 


Integer  array  dimensioned  to  16  and  containing  the 
unpacked  y-axis  titles.  The  titles  are  unpacked  from 
four  to  a  word  (in  array  LABY)  to  one  to  a  word 
( LABY2)  in  ISN  72  of  entry  point  INITP  (subroutine 
FHPLOT).  Unpacking  is  performed  only  after 

translation  from  EBCDIC  to  ASCII  has  taken  place. 


Integer  array  dimensiond  to  five  vcrds,  containing 
the  number  of  characters  in  each  of  the  five  x-axis 
titles.  the  index  (i.e.,  subscript)  of  LEN  required 
to  retrieve  the  length  of  a  specific  title  is  the 
value  of  the  variable  IXAX  (see  COMMON  block  XAXPAR) 
for  that  title.  LEN  is  assigned  values  in  the  BLOCK 
DATA  subroutine. 

Integer  variable  whose  value  is  the  number  of 
characters  in  the  y-axis  title. 
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APPENDIX  D. — LIBRARY  ROUTINES 

All  library  subroutines  and  functions  required  by  the  program  are 
described  in  the  following  appendix. 


GENERAL-PURPOSE  ROUTINES 


ELT 1 


Initializes  CPU  timer  to  0  for  determination 
of  CPU  time  used  in  program  execution. 


ELT3 
( ITIME) 


Places  CPU  time  used  since  last  call  to  ELT  1 
into  12-byte  array  ITIME  in  format 
HH :MM : SS. FFF. 


TREAD  ( INPUT, 
INLEN,  LENBUF) 


PGMASK 


Reads  characters  from  terminal  into  array 
INPUT.  INLEN  is  the  number  of  characters 
read,  and  LENBUF  is  the  length  of  INPUT  in 
bytes.  If  the  number  of  characters  typed  is 
greater  than  LENBUF,  then  the  input  will  be 
truncated  to  length  LENBUF,  and  an  error 
message  will  be  printed. 

Eliminates  error  message  when  exponent 
underflow  occurs. 


TREA  ( LEN , IBUF ) 


UNPAK  (LEN,  IBUF 1 
IBUF2 ) 


Translates  LEN  characters  in  array  IBUF  from 
EBCDIC  to  ASCII  character  codes.  Characters 
must  be  packed  one  to  a  byte  in  IBUF. 

Unpacks  LEN  characters  in  array  IBUF 1  and 
places  unpacked  characters  in  array  IBUF2  for 
use  by  plotting  routines.  Characters  must  be 
packed  one  to  a  byte  in  IBUF  1 ,  and  will  be 
placed  in  IBUF2  one  to  a  word,  right- 
justified.  IBUF2  must  be  dimensioned  to  at 
least  LEN  words. 


TWRITE  (IBUF,  LEN)  Writes  LEN  characters  in  array  IBUF  to 

terminal,  without  preceding  or  succeeding 
carriage  return/line  feed.  Characters  in  IBUF 
are  packed  one  to  a  byte.  This  routine  is 
used  mainly  for  input  prompting.  Example: 
CALL  TWRIIE  ('INPUR  LAMBDA:'  14) 

IREAD  (INTEG)  Reads  a  single-integer  number  from  the 

terminal.  This  routine  uses  TREAD  to  read  the 
terminal  input,  and  then  uses  FORTRAN  integer 
conversion  routines.  The  integer  input  is 
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TREADC  ( INPUT, 
INLEN,  LENSUF) 

FREAD  ( FLOAT) 


Q  ( A,  B) 

rO( ARG,  RESULT) 

NL9  (SUBR,  XL, 

XH,  ACCUR,  INT, 
LIM,  RESULT, 
ERROR,  NFUN ,  IER ) 

LOGCOM  (I,  J) 

COMB  (I,J) 

BINOM  (I,J,P) 

ERFC  (X) 


placed  in  the  variable  INTEG  on  return.  If 
the  character  '#'  is  entered  as  the  first 
character  from  the  terminal ,  IREAD  executes  a 
STOP  command,  halting  execution  of  the  program 
and  returning  control  to  the  operating  system. 

Same  as  TREAD,  except  all  lower-case 
characters  in  INPUT  are  forced  to  upper  case. 

Reads  a  single  floating-point  number  from  the 
terminal.  This  routine  uses  TREAD  to  read  the 
terminal  input,  and  then  uses  FORTRAN 
floating-point  conversion  routines.  A  decimal 
point  or  "E"-type  exponent  is  optional.  The 
floating-point  number  is  placed  in  the 
variable  FLOAT  on  return.  If  the  character 
"#"  is  entered  as  the  first  character  from  the 
terminal,  FREAD  executes  a  STOP  command, 
halting  execution  of  the  program  and  returning 
control  to  the  operating  system. 

Function  which  returns  the  value  of  Marcum's 
Q-function  for  arguments  A  and  B. 

Computes  the  zero-order  modified  Bessel 
function  of  ARG,  returning  the  result  in 
RESULT. 

Estimates  the  integral  of  a  function 
(calculated  in  subroutine  SUBR,  supplied  by 
user)  between  limits  XL  and  XH,  and  returns 
estimate  in  RESULT. 

Function  which  computes  logg  (I,J)  where  (I,J) 
=  I ! / [ J  !  ( X  — J ) ! ]  .  The  real  *4  result  is 
returned  as  the  value  of  the  function 
LOGCOM.  LOGCOM  must  be  declared  REAL  *4. 

Function  which  computes  (I,J)  =  I!/[J!(I- 

J)!].  Ibe  real  *4  result  is  returned  as  the 
value  of  the  function. 

Function  which  computes  the  binomial 
coefficient  BINOM  =  ( I , J) PJ  (1  -  P)I-J. 

FORTRAN  library  function  which  computes  the 
complementary  error  function  for  argument  X. 
Idle  result  is  returned  as  the  value  of  the 
function . 
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FACT  (I) 

PLOTTING  ROUTINES 
INITT  (IBAUD) 

ERASE 

ANMODE 

MOVABS  (  IX , IY ) 

B INITT 

XFRM(I) 

YFRM(I) 

XMFRM(I) 

YMFRM(I) 

YTYPE(I) 

NPTS ( INPT) 


Function  which  return  FACT(I)  =  I!. 

For  Tektronix  graphics  terminals  only. 

Initializes  PLOT-10  library.  IBAUD  is  the 
character  transmission  rate  at  which  the 
terminal  will  communicate  with  the  host 
computer  (in  characters  per  second).  Must  be 
called  exactly  once  during  program  execution. 

Erases  the  screen  of  the  Tektronix  terminal. 

Converts  the  terminal  to  the  alphanumeric 
mode . 

Moves  the  cursor  in  visibly  to  the  absolute 
screen  position  represented  by  the  coordinates 
(IX,IY).  This  routine  is  used  in  order  to 
position  the  cursor  at  the  desired  starting 
point  of  a  visible  line  or  character  string. 

Resets  PLOT-10  option  flags  back  to  their 
default  values. 

Selects  options  for  x-axis  major  grid  lines 

and  tic  marks  based  on  value  of  argument  I. 

Selects  options  for  y-axis  major  grid  lines 

and  tic  marks  based  on  value  argument  I. 

Selects  options  for  x-axis  minor  tic  marks 

based  on  value  argument  I. 

Selects  options  for  y-axis  minor  tic  marks 

based  on  value  of  argument  I. 

Selects  y-axis  scale  type  (i.e.,  linear, 
logarithmic,  etc.)  based  on  value  of  argument 
I.  For  1  =  2,  logarithmic  scaling  is 

selected.  Default  is  linear  scaling. 

Passes  to  the  plot  library  the  number  of  (x,y) 
data  points  to  be  plotted  on  a  given  plot. 
INPT  is  the  number  of  points  that  are  passed. 
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LINE ( I ) 


DLIMX(XMIN, XMAX) 


DLIMY(YMIN, YMAX) 


CHECK  ( X,  Y) 

DISPLAY ( X , Y ) 

CPLOT ( X , Y ) 

JUSTER  (LEN, 
LABEL,  OUST, 
IFILL,  LENF, 
IOFF) 


NOTATE  (IX,  IY, 
LEN,  LABEL) 


Selects  line  type  to  be  used  by  the  plot 
library  when  connect ing  the  points  being 
plotted.  1=0  selects  a  solid  line;  1=2 
selects  a  dashed  line. 

Passes  to  the  plot  library  the  minimum  and 
maximum  values  of  X  coordinates  of  the  points 
to  be  plotted. 

Passes  to  the  plot  library  the  minimum  and 
maximum  values  of  the  Y  coordinates  of  the 
points  to  be  plotted. 

Scales  the  data  to  fit  in  the  plotting  window, 
calculates  x-  and  y-axis  labels  and 
increments,  performs  other  checking 

functions.  X  and  Y  are  the  arrays  containing 
the  X  and  Y  coordinates  of  the  data  to  be 

plotted.  The  number  of  points  in  X  and  Y  to 
be  plotted  are  specified  by  calling 
subroutine  NPTS. 

Plots  the  (x,y)  data  on  the  screen,  and  draws 
the  axes  and  axis  labels. 

Plots  an  additional  set  of  (x,y)  data  points 
on  an  existing  plot,  scaling  the  new  data  to 
the  existing  plot  and  clipping  points,  if 
necessary.  X  and  Y  are  arrays  containing  the 
X  and  Y  coordinates  of  the  additional  data 

points . 

Left-,  right-  or  center- just i f ies  (selected 
by  JUST)  the  character  string  of  length  LEN  in 
array  LABEL  (packed  one  character  to  a 
word).  The  length  of  the  character  string 

without  fill  characters  (LENF)  and  the 
distance  in  screen  raster  units  between  the 
justification  point  (e.g.,  center)  of  the 

string  and  the  starting  point  of  the  string 
(IOFF)  are  returned.  This  routine  is  used  to 
center  x-axis  titles. 

Writes  the  character  string  of  length  LEN  in 
array  LABEL  ( packed  one  to  a  word)  at  screen 
location  (IX,IY).  This  routine  is  used  to 
write  x-axis  titles. 
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VLABEL  (LEN, LABEL) 


IBASEX ( I ) 

IBASEY ( I ) 
COMGET  ( I BASE) 


>*■ 


Writes  the  character  string  of  length  LEN  in 
array  LABEL  (packed  one  to  a  word)  vertically 
(i.e.,  a  line  feed  and  backspace  are  performed 
between  each  character).  The  writing  of  the 
vertical  label  begins  at  the  current  cursor 
position,  and  may  be  specified  by  performing  a 
MOVABS  immediately  prior  to  the  call  to 
VLABEL.  This  routine  is  used  to  write  y-axis 
titles . 

Function  which  return  as  its  value  the  address 
of  the  PLOT-10  x-axis  COMMON  variable 
specified  by  the  value  of  the  argument  I. 
Used  to  retrieve  previously  set  PLOT-10 
options  and  variable  values. 

Same  as  IBASEX,  except  for  y-axis  variables. 

Function  which  returns  as  its  value  the  value 
of  the  COMMON  variable  whose  address  is  the 
value  of  the  argument  IBASE  (obtained  by  call 
to  IBASEX  or  IBASEY).  Used  to  retrieve  PLOT- 
10  values  such  as  axis  lengths  and  locations. 
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APPENDIX  E.— VARIABLE  DEFAULT  VALUES 

The  following  appendix  names ,  describes,  and  gives  the 
values  for  certain  program  variables. 

default 

Variable 

Name 

Variable  Description 

Default  Value 

A 

Frequency  hopping  channel  overlap 
factor. 

1. 

B 

Signal- jamming  center  frequency 
difference  factor. 

2. 

BETA  1 

Signal-to-noise  ratio  per  word  in  decibels 

.  13. 

BW 

Frequency  hopping  bandwidth  in  Megahertz.- 

25. 

D 

PN  receiver  loss  factor  (also  called 
the  receiver  time-bandwidth  product) . 

1. 

FB 

Information  bit  bandwidth  in  Kilohertz. 

25. 

FCM 

Chip  (coded  bit)  bandwidth  in  Kilohertz. 

25. 

G 

PN  processing  gain  in  decibels. 

30. 

GAMMA 

Jamming- to- signal  ratio  in  decibels. 

0. 

NBPC 

Chips  (coded  bits)  per  word 

1 

NBPW 

Word  length  in  bits 

1 

NERR 

Chip  error  threshold  (the  number  of  bit 
errors  required  in  a  coded  word  in  order 
to  cause  a  word  error). 

1 

N  JAM 

Number  of  jammed  channels. 

1 

UPBJ 

Portion  of  band  jammed  (partial  band  jamming).  0.25 

t /-.i.  ^i,»u*R»«Cr 
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